Quantum tunneling of superconducting string currents. 
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Abstract 

We investigate the decay of current on a superconducting cosmic string 
through quantum tunneling. We construct the instanton describing tunneling 
in a simple bosonic string model, and estimate the decay rate. The tunnel- 
ing rate vanishes in the limit of a chiral current. This conclusion, which is 
supported by a symmetry argument, is expected to apply in general. It has 

important implications for the stability of chiral vortons. 
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I. INTRODUCTION 



Our universe may contain cosmic strings which are relics of a symmetry breaking phase 
transition at early times. (For a review, see Q.) In many models these strings can be 
superconducting ||, and the resulting current can have important effects on the string 
dynamics. In particular, if the current is large enough, the pressure due to the charge 
carriers can oppose the tension of the string and lead to a stable "vorton" state ||. 

The string current depends on the winding number density of the phase of a condensate 
field; in a loop the total winding number will be an integer, and thus classically cannot 
change as long as this phase is well-defined. When the current becomes large enough, there 
is the possibility of a "quench", where the condensate is driven to 0, and thus the phase can 
unwind and current can be lost. However, even below this threshold current, there is the 
possibility that current can be lost through quantum tunneling effects. It is those effects 
which we consider here. 

The simplest field theory model that exhibits bosonic superconductivity is of the form, 

C = D^D^a + d»<j> - -J^r - V(a, 0), (1) 

where 

f iiv — d fJi Z u — d u Z fJi , 
D„a = (fy - xeZ^a . (2) 

This model has a local U(l)i oca i symmetry associated with the gauge field Z M and a global 
symmetry U(l) g i b a i of the phase transformations of the field 0. In order to have a super- 
conducting string model, we should also impose that the interaction potential between the 
two complex scalar fields, <fi and a is such that the broken U(l)i oca i symmetry gives rise 
to a cosmic string configuration and the U(l) g i i, a i becomes broken only inside the string. 
A bosonic condensate will then be formed on the string, with massless Goldstone bosons 
playing the role of "charge carriers" . Note however that in this model the condensate field 
cf) is not coupled to any gauge field, so the current is neutral. 
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Low energy excitations of the condensate along a straight string in the z direction will 
have the form 

<f>(x, y, z, t) = <p(x, y) exp[i9(t, z)\ . (3) 

Using this ansatz, we can write an effective action for the field 9 by inserting Eq. (0) in the 
four dimensional action to get, 

s = s J dzdt[{d t ef - {d z ef} (4) 

where 

£ = J dxdy\(j)\ 2 . (5) 

We can now define the neutral current associated with the field living on the string 
worldsheet by[] 

J a = d a 9, (6) 

so the current is a measure of the winding number per unit length of the phase of the bosonic 
condensate on the string. 

If we imagine the straight string as part of a large loop, we can write the total winding 
number 

N = ^-<fd a 9da, (7) 

Z7T J 

which is classically conserved. Nevertheless, this is not the whole story. It was already 
noticed by Witten in his original paper on superconducting strings that quantum effects 
would be able to decrease the current. He argued that since the field would have to unwind 
in order to decrease the value of N, the process would have to involve some kind of localized 
region on the string worldsheet where the absolute value of the field <fi would vanish. 



1 Here and below, Latin indices take values and 1, Greek indices take values from to 3, and 
the signature of the metric is goo > both in four dimensions and on the string worldsheet. 
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Several authors have estimated the decay rate of the current following similar arguments 
In this paper we shall study the quantum tunneling of the current by analytically 
computing the worldsheet instanton for this process and the corresponding decay rate in a 
loop with a uniform current. 



II. WORLDSHEET INSTANTON 

As a first approximation we will study the process on a straight string and will neglect 
any backreaction effect on the profile of the string. Therefore, we will only be concerned with 
the evolution of the field in a fixed background. Also, we will not consider situations where 
the unwinding of the field takes place differently at different locations in the plane transverse 
to the direction of the string. This enables us to integrate over transverse directions, and 
use only the effective Euclidean action for the field 9, 

S E = T,Jdz dr[(d T 9) 2 + (d z 9) 2 ] (8) 

The Euclidean field equation for 9 is then 

d a d a 6 = 0. (9) 

The limits of validity of this approximation will be indicated below. 
Let us first assume that the current is spacelike, 

j a r < o. (io) 

Then we can choose a frame of reference where the charge density vanishes, Jq = (pure 
current). We are interested in studying the quantum tunneling process from an initial state 
of homogeneous pure current of winding N, 

* = ^, (id 

where L is the length of the loop, to a final state of winding N — 1. The corresponding 
initial current is therefore, 



Semi classically, the tunneling is described by a Euclidean instanton solution 9{r, z) which 
approaches the initial configuration ( |TTD at r — > ±00. The field configuration immediately 
after tunneling is given by 6(0, z). Since the winding along the z-axis r = should be one 
unit less than that at r — > ±00, it is clear that the field 6(r, z) should have a vortex below 
and an anti- vortex above the 0-axis (or vice versa). The corresponding solution of Eq. @ 
is 

/./ ( * \ ( z \ 2llN / \ 

0(T,z) = arctan — arctan H — z, (13) 

v ' \T-aJ \T + aJ L v ' 

where 2a is the vortex separation. This solution is well known from hydrodynamics , where 
d a 8 corresponds to the velocity of an incompressible fluid. The solution describes a vortex- 
antivortex pair on the background of a fluid of constant velocity. The vortex separation a is 
chosen so that the two vortices are in equilibrium due to the balance between their mutual 
attraction and the oppositely directed Magnus forces exerted by the background. This gives 

1 

a = —. (14) 

The equilibrium is unstable, as it should be, since the instanton solution should have a 
negative mode . 

For our purposes, instead of the variable 6 it will be convenient to introduce the dual 
field x defined by 

d a X = £abd b 6. (15) 

Like 6, x obeys Laplace's equation, d a d a x — away from the vortex cores, but a vortex in 
9 becomes a source for x- This can be used to write the equation for x m the whole space, 

d a daX = 27rp, (16) 



where the "charge density" p is given by 

p = 5(x-x + )-5(x-x_). (17) 

Here, x is a vector in the (r, z)-plane, x± = (±a, 0) are the locations of the vortex centers, 
and S(x) is a delta-like function normalized so that 

5{x)d 2 x = 1. (18) 

The width of the function 5(x) is comparable to that of the vortex core and its precise shape 
will not be important for our considerations. We thus have an analogy with 2-dimensional 
electrostatics. The problem to be solved is that of two opposite charges of magnitude 1 in 
a constant external field of magnitude 

E = (2-rrN/L, 0). (19) 

The corresponding solution of flUf) outside the sources is 



x = In X X+ + E • x. (20) 



x — x 



The mutual attraction of the charges must be balanced by the force of the external field 
acting on them in opposite directions, which occurs when E = I /2a, in agreement with (|i~i|). 
See Fig |. 

The Euclidean action corresponding to the field equation (|16D can be written as 

S E = -Y,j d 2 x[(d aX ) 2 + ^p X }. (21) 

The overall factor has been chosen so that the actions @ and (|21|) coincide when x is a 
solution of (|T6| ) and 6 is related to x by Eq. (|T5|). 
The current decay rate is given by 

T = Ae~ B , (22) 

where 



FIG. 1. This figure shows the lines of constant x f° r the Euclidean instanton solution Eq. 
( |l3| ) on the string worldsheet. The attractive force exerted between the vortex and anti-vortex is 
balanced by the Magnus force due to the background current. 

B = S E -SP, (23) 

Se is the instanton action and is the action for the constant current configuration flnp. 
The semiclassical approximation applies when B ^> 1. In this regime, the tunneling rate is 
determined mainly by B, and the value of the prefactor A is not very important. 

A. Evolution subsequent to the decay. 

It is easy to check that our solution (|i~3"l) conserves energy and momentum as well as 
charge on the worldsheet. We can also analytically continue it into Lorentzian time to 
obtain the evolution of the current after the tunneling process has occurred. By calculating 
the charge density on the Lorentzian worldsheet after the current decay, we can see (Fig. 
0) that this solution represents the propagation of two lightlike trains of opposite charge 
travelling in opposite directions along the string. They leave behind a region of the same 
winding number density as the initial conditions, but if we now calculate the total winding 
along the string we see that 



Winding density 
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FIG. 2. Evolution of the winding and charge density in the string condensate after the quantum 
current decay process. The two travelling waves of opposite charge fly out from the center of the 
nucleation process and their negative winding with respect to the positive background decreases 
the overall integrated winding along the string. 



This is because the two travelling waves have each negative winding number density, in 
other words they have the opposite winding to the original current. 



In order to calculate the decay rate for this process we have to compute the Euclidean 
action for the instanton, or in terms of the electrostatic analogy, the difference in energy 
between the pure background configuration and the equilibrium solution described above. 
However, both energies are infinite, due to the infinite two-dimensional volume occupied by 
the background electric field. Moreover, one finds that a formal large-distance cutoff gives 
different answers, depending on the details of the cutoff procedure. We have therefore intro- 
duced a "physical" cutoff by replacing the homogeneous electric field with a field produced 
by a finite distribution of charges, p exi (x). The solution (^) is then replaced by 




(24) 



III. DECAY RATE 



S 



X = In — ^7 + Xext{*) = X( x ) + Xext(x), (25) 
|x — x_| 

where Xext( x ) is the solution of ([R]) with the source p ex t(x) and it is assumed that Xext( x ) 
vanishes at |x| — > oo. It is assumed also that the external charges are located at distances 
|x| ^> a, and in the end we take the limit when these charges are removed to infinity, while 
keeping the electric field near the origin fixed and equal to (JT9|) , 

- J d 2 xp ext {*)^ = E. (26) 

We also have to deal with the short- distance logarithmic divergence in the self-energy of 
the two charges. This divergence is an artifact of our approximation treating the vortices 
as pointlike objects. It is easily fixed by replacing InO by In 5, where 5 is the characteristic 
radius of the vortex core. Although 5 is defined only by the order of magnitude, the result 
is not sensitive to the precise definition, due to the logarithmic character of the divergence. 

After integration by parts in Eq. (pT|), the energy difference between the instanton con- 
figuration and the background can be expressed as 

B = -2ttS (J pxd 2 x - J pextXextd 2 x^j = -2ttS[x(x + ) - x(x_)] - 2vrS J p ext (x)x(x). (27) 

Since the support of p e xti^) is at |x| ^> a, we can use the large- distance expansion of x, 

2x • a , . 

(28) 



in the last term of ( p27D . With the aid of Eq. (pq) , this gives 

B = 47r£[ln(2a/<f) -2Ea), (29) 

where 5 is the vortex core radius, as explained above. Note that the result ( p^| ) is insensitive 
to the details of the external charge distribution p ext (x). 

The action ( p5| ) is maximized for ao = 1/2.E, in agreement with (p!4|). It is easily verified 
that d 2 S/da 2 (ao) < 0, indicating that a = ao is a local maximum of the action. In terms of 
the electrostatic analogy, this indicates that a = ao is an unstable equilibrium, as it should 



be for an instanton. Substituting a = ao in (|2Q ) we obtain finally 



B = 47rE[ln(L/27riV5) - 1], (30) 

and for the tunneling rate 
where K = 47rS. 

The prefactor A for tunneling in (1 + 1) dimensions is given by A = BD, where D is 
a determinantal factor defined in Ref. @ and the factor B comes from including the zero 
modes contribution to the determinant. 

The exact calculation of this prefactor is a complicated task and we do not attempt it 
here. Our main interest here is in the dependence of the prefactor on the worlsheet current. 
In the Appendix we show that we should expect this prefactor to scale as, 



A~ B5~ 



2 ty 2 - (32) 



\2tiNSJ 

Analytic estimates and numerical calculations show || that, if the fields a and <fi are 
weakly coupled to each other and to themselves, then E > 1 in a large portion of the 
parameter space of the model (|I|). In this case, K ^> 10, and it follows from Eqs. ( |3"T| ) and 
( p2f) that the current decay is strongly suppressed, as long as the condensate winding does 
not get very large, namely, iV <C L/27r5. 

We finally discuss the validity of the approximations we made in deriving the tunneling 



rate (31). First, we assumed that the vortex cores give a negligible contribution to the 



instanton action, compared to the contribution coming from the region 

5 |x — x±| < a. (33) 

This is justified if the logarithm in Eq. fl29|) is large, 

\n(L/2TtN5) > 1, (34) 

that is, when the current is very small. We expect the vortex radius 5 to be comparable to 
the thickness of the condensate; then Eq. ([34]) says that the wavelength of the condensate 
should be large compared to its thickness. 
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We also assumed that the instanton solution can be factorized as in (|3]), with 9 inde- 
pendent of the coordinates (x, y) in the plane transverse to the string. This approximation 
is likely to break down in the vicinity of the vortex cores, where 6 varies on a length scale 
comparable to the string thickness, but we expect it to be valid in the range (|33| ) which, 
assuming fl3"3|), gives the dominant contribution to the action. Eq. (Q) is thus the only 
condition we need to impose to justify the approximations made above. 

IV. CHARGED CURRENT INSTANTON. 

We would like now to extend the previous analysis to the case where the condensate is 
coupled to the electromagnetic field A^. The simplest way to achieve this is to gauge the 
U (1) global symmetry we have been discussing. In this case the action for the model takes the 
form, 

S el = J ' d A x^D»a --J^r + V^V a (j>-V{aA)-\F llv F^ (35) 

where 

F^v = d^A u — dyA^ , 

= (fy - icA^cj) . (36) 

As in the previous case, we are interested in the quantum mechanical decay of a spacelike 
current of the form 

<p = <f>(x,y) exp[z6»(r, z)\ . (37) 

In terms of the angular variable the effective action becomes 

S el = S cff J dzdr(d a 6 + eA a f - 1 J >x . (38) 

where S e fj can be written in terms of the parameters in the original U(l)i oca i x U(l)i oca i 
Lagrangian. This expression is correct if the scalar field condensate and the vector field 
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remain approximately constant in the transverse direction to the string. This is ensured if 
the penetration depth of the vector field in the core of the string is large compared to the 
condensate thickness 

The important difference in this case is that the electromagnetic field is not confined to 
the worldsheet, so our instanton is somewhat more complicated. Adding the electromagnetic 
field makes the Euclidean equations of motion for the worldsheet look like those of a two 
dimensional superconducting sheet embedded in a four dimensional Euclidean space. Our 
initial condition at r = — oo is a sheet carrying a constant planar current with its accompa- 
nying magnetic field. The instanton includes a vortex-antivortex pair in equilibrium under 
the action of the Magnus force, the electromagnetic force due to the long-range magnetic 
field of the vortices, and (to a smaller extent) the force due to the scalar field distribution 
on the sheet. 

When this work was in progress we learned that this type of instanton has been discussed 
by Duan |9j in the context of a real superconducting wire. The instanton solution cannot 
be found analytically in this case, but remarkably the instanton action, (excluding the 
contribution from the vortex cores) can still be found exactly M, 



B 



4vrS 



a 



In 



1 + a In 



2a 

T 



/47riVa\ 



(39) 



where a = — , and N, L, and S have the same meaning as before. Note that this reduces 
to our Eq. (|29|) in the neutral limit e — > 0. The action is maximized at 



a 



(—)- 

\47rNJ i 



1 



+ a In 



2ttN8 



(40) 



and finally it becomes 



B 



4ttS 



a 



In 



a In 



L 



2nN5 



a 



1 + a In 



2ttN8 



(41) 



As we discussed at the end of Section III, we are interested in the range of parameters 
where S > 100 and \n(L/2nN5) > 1. With e 2 ~ 10~ 2 this gives a > 1 and Eq. ([|]]) reduces 
to 
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Following the discussion for the neutral case we see that the tunneling rate for the charged 
case is 

r ~ A { a H^m)) (43) 

which clearly indicates that this process is very much supressed in the electromagnetic case 
as well. 



V. DISCUSSION 

Our result ([31]) was derived in the frame where the global charge density vanishes, Jo — 0. 
A generalization for an arbitrary frame is easily obtained if we notice that the rate T is 
Lorentz invariant, so we only need to express it in terms of the invariant quantity J a J a - 
With the aid of Eqs. (|3~1"D, fl3"2|) and assuming (|34]), we obtain 

T ~ K5- 2 X^ K - 2 ^ 2 ln(l /X), (44) 

where 

X=\J a J a 5 2 \. (45) 



It is interesting to note that the rate fl44f) vanishes in the chiral limit J a J a = 0. This is 
not surprising: since T is an invariant, it can only depend on the invariant combination J a J a 
and should therefore vanish for a chiral current. We expect this conclusion to be valid for 
strings with charged as well as neutral currents. The stability of chiral vortons with respect 
to this type of decay reinforces the argument that cosmic string theories able to produce 
chiral vortons will have serious constraints from their overproduction in the course of the 
history of the string network |T0|-[T^ . 



For a timelike current, J a J a > 0, instantons of the type we discussed here do not exist. In 
this case, one can go to the Lorentz frame where the current vanishes and there is only a non- 
zero global charge density. Charge is conserved on the worldsheet and cannot therefore be 
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reduced by worldsheet instantons or any other processes operating on the string worldsheet. 
It is possible that when the charge density gets high enough, charges can escape from the 
string by tunneling into the external space. We shall not consider such tuneling processes 
in the present paper. 



We would like to thank Jaume Garriga for very helpful conversations. This work was 
supported in part by funding provided by the National Science Foundation. J.J.B-P was 
supported in part by the Relativity Group PPARC Rolling Grant. 



In this Appendix we will estimate the dependence of the decay rate prefactor on the 
uniform current on the worldsheet. The usual procedure for this type of calculation involves 
the computation of the imaginary part of the energy coming from the contribution around 
the one-instanton solution of the Euclidean equations of motion. The decay rate per unit 
volume per unit time is given by @, 



where S E and S" E denote the second order perturbation expansion of the instanton and the 
background Euclidean action respectively. In our notation we have 



Furthermore, the first term, proportional to N^, is included to take into account the 
normalization of the zero modes, and in general it can be shown to be proportional to the 
instanton action B [[5J. Consequently we do not have to consider those in the following 
calculation. This is signaled by the expression Det' which denotes the determinant of the 
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(46) 



q q(0) _ O 

D E — O e — D . 



(47) 
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fluctuations once the zero modes have been excluded. In our system the zero modes corre- 
spond to the collective motion of the vortex pair in the (r, z) plane. It is clear that translating 
them together does not alter the action, so this kind of fluctuations must correspond to zero 
modes. 

Among the modes included in Dei! is the negative eigenvalue mode of the fluctuations 
around the instanton solution. As we discussed earlier, this negative mode corresponds to 
the fluctuation that changes the distance between the vortices in the direction perpendicular 
to the current flow. 

There is also another type of fluctuation which corresponds to a relative motion of the 
vortices in the direction parallel to the current. This is a positive mode. 

If the current is very small, in other words, if the distance between the two vortices 
is very large compared to the size of the vortex core, it could be expected that the only 
fluctuations which would not cancel with the ones in the background are the fluctuations 
associated with motions of the vortices. In this case the effective one-loop contribution to 
the prefactor can be written as the following integral, 



where the factor of B is the usual normalization factor for the two zero modes in 1+1 
dimensions. This expression singles out the contribution from the modes that are related to 
the position of the vortices by defining z r and r r as the relative distance between the vortices 
on the worlsheet coordinates. Furthermore, ao denotes the distance along the Euclidean time 
direction that maximizes the instanton action, (5 = 47rS and /i is a normalization factor of 
order |0| ~ S^ 1 summarizing the result of doing the small fluctuation integral around the 
one vortex configuration. 

Performing the gaussian integrals we arrive at 




(48) 




(49) 
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